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ABSTRACT 

We present the evolution of the radio emission fro m the 2.8-s pulsar of th e double pulsar system 
PSR J0737— 3039A/B. We provide an update on the Burgav et al. ( 2005f ) analysis by describing 
the changes in the pulse profile and flux density over five years of observations, culminating in 
the B pulsar's radio disappearance in 2008 March. Over this time, the flux density decreases 
by 0.f77 mJy yr _1 at the brightest orbital phases and the pulse profile evolves from a single to 
a double peak, with a separation rate of 2.6° yr _1 . The pulse profile changes are most likely 
caused by relativistic sp in precession, but can not be easily explained with a circular hollow-cone 
beam as in the model of Clifton fc Weisberg] (2008). Relativistic spin precession, coupled with an 



elliptical beam, can model the pulse profile evolution well and the reappearance is expected to 
happen in ~ 2035 with the same part of the beam or in ~ 201 4 if we assume a symmetric beam 
shape. This particular beam shape predicts geometrical pa rameters for the tw o bright orbital 
phases which are consistent and similar to those derived by iBreton et al. I (l2008h . However, the 



observed decrease in flux over time and B's eventual disappearance cannot be easily explained 
by the model and may be due to the changing influence of A on B. 
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Introduction 

Double neutron stars are systems of two neu- 
tron stars orbiting one another. Only ten such 



syste ms are known so far ([Stairs 2008; K eith et al 
2009). The detection of both neutron stars as ra- 



dio pulsars is a rare situation, as it requires the 
electromagnetic beam from both pulsars to in- 
tercept our line-of-sight during the short interval 
when both pulsars are active. The only such sys- 
tem, J0737— 3039A/B, was discovered in late 2003 
with the 64-m Parkes radio telescope in Australia 
as part of a high -latitude multibeam survey of 
the S outhern sky (jBurgav et al. 1 120031: iLvne et all 
2004). It provides an amazing laboratory for the 
study of relativistic gravity and the most pre- 
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cise test to date of general relativity (G R) in the 
strong- field regime ( Kramer et al. 2006| ) . The two 
pulsars orbit each other in a 2.4-hr orbit, the short- 
est of any observed double neutron star binaries, 
with moderate orbital eccentricity of 0.088. 

The first-born pulsar is believed to have been 
recycled, hereafter referred to as 'A', has a short 
spin period of 23 ms and the second-born com- 
panion, hereafter 'B', spins with a longer 2.8 s 
period. With their strong gravitational fields 
and rapid motions, the system is expected to 
show large relativistic effects and, indeed, has 
a relativistic advan ce of periastron of 17° yr _1 
( Kramer et al.ll2006l ). the highest rate of any ob- 
served pulsar binary. Another of these effects 
is the relativistic spin precession, also known as 
geodetic precession, of the pulsar's spin axis about 
the orbital angular moment um. The theoreti- 
cal eq uatio n formulated b y Barker fc O'Connelll 
(|l975l ) and iBoerner et ail (|l975l) predicts a spin 
precession rate of 5.061(2)° yr" 1 for pulsar E0. 
Recently this rate has been constrained to be 
4.8(7)° yr _1 by fitting a model to the chang- 
ing morphology of the eclipses of A seen when 
it passes behind the a bsorbing magnetosphere of 
B ([Breton et al.l 12008). The spin precession rate 
predicted by general relativity means that it takes 
75 yrs for B's spin axis to complete a full preces- 
sion cycle about the orbital angular momentum. 
This high precession rate is expected to have dra- 
matic effects on B's pulse properties as a result 
of the change in apparent geometry with respect 
to our line-of-sight. These include the evolution 
of its pulse profile, l ight curves, and flux densities 
(jBurgav et al.ll2005h . 

The orbital plane of this system is almost edge- 
on to us. The orbital inclination was calculated to 



be 88.7° through timing ([Kramer et al.ll2006j ) and 
89.7° through modeling of the interst ellar scintil 



lation properties of the two pulsars (|Coles et al 



2005). Eclipses of A by B are a fortunate conse- 



quence of this favorable geometry. The 2.9 lt-s di- 
ameter of the orbit, determined through radio tim- 
ing, and the 30-s duration of the eclipse, show that 
the magnetosph ere of B is roughly 10% of its light- 
cylinder radius (|Lvne et al.1120041 ). The pulse pro- 



file of the more energetic pulsar A is stable. The 



limit on variations of A's profil e width at 10% of 
the p ulse peak is only 0.5° yr _1 ([Manchester et al 



2005). This can b e explained by a range of g eomet 



rical parameters. [Manchester et al. ( 20051) found 



that both circular and elliptical beams could fit 
the da ta, with a slight prefe rence for an elliptical 
beam. Ferdman et al. (120081) found that the mis- 



alignment between spin and orbital angular mo- 
mentum axes of pulsar A is very low, with a 95% 
upper limit of 14°, assuming emission from both 
magnetic poles. This supports the idea that the 
formation process of pulsar B was relatively sym- 
metric. 

The high spin-down energy loss rate of A likely 
causes the magnetosphere of B to be deformed in 
a way similar to that of the Earth by the Solar 
wind ([Lvutikovl 120041 ). The long tail will point 
towards the Earth around the eclipse of A by B, 
near the inferior conjunction of B (i.e. at an or- 
bital phase of 270°, as measured from the ascend- 
ing node) . We find that the bright emission from B 
is detected in two different orbital phase windows: 
at -180°-240° (hereafter BP1), and ~260°-310° 
(hereafter BP2). In addition to these two bright 
phases, the pulsar o ften shows weak emission at 
other orbital phases ([Lvne et alj|2004l) . The exact 
explanation of this exc eptional behavio r is chal- 
lenging. The model of iLvutikovl (|2005[ ) explains 
the phenomenon generally, but fails to produce the 
second bright phase at the correct orbital phase 
(~245°-270° instead of -260°-310°). It also can 
not explain the evolution of the two bright phases 
(see §2.3). 

Another effect of A's electromagnetic radiation 
on B's magnetosphere is drifting subpulses in the 
emission of B which appear si milar to those seen 



for so me normal radio pulsars ( McLaughlin et al 



2004). The drifting feature can be seen only in 



1 Here, and throughout the paper, the number in parentheses 
is the 1-a uncertainty in the last quoted digit. 



BP1 and not in BP2. The separation of the drift 
bands is equal to the period of A, showing that the 
drifting features are a direct result of the impact 
of A's 44 Hz electromagnetic radiation on B. A 
possible explanation for the occurrence of drifting 
features only in BP1 is the particular geometric 
configuration of t he two pulsars wit h respect to 
our line-of-sight. iFreire et al. (2009) proposed a 
new technique for timing the double pulsar system 
by describing a geometrical model which predicts 
the delay of B's subpulses relative to A's radio 
pulses. 
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In this paper, we present a study of the extreme 
variation of B's pulse profile over time, which con- 
cludes with its total radio disappearance in March 
2008. We present an analysis based on more than 
five years of data taken with the Green Bank Tele- 
scope. In Sj2]we present our observational data and 
our analysis methods. We also discuss the disap- 
pearance of B's emission and the evolution of the 
emission from the beam. In |j3l we estimate the ge- 
ometrical parameters of pulsa r B by using the em- 
pirica l circular beam model of Clifton fc Weisberg] 
(|20081 ) for the two-dimensional pulse profile of a 
precessing pulsar. Then we extend this model to 
an elliptical beam that explains the profile evolu- 
tion well and results in realistic geometrical pa- 
rameters. In [21 we discusses how the different 
types of beam shapes relate to the observed pul- 
sar data and show that this is the best explanation 
for the behavior of B. 

2. Observations and Analysis 

We observed J0737-3039B with the 100-m 
Green Bank Telescope (GBT) in West Virginia 
from 2003 December 24 to 2009 June 20. Our ob- 
servations are at multiple frequencies - because 820 
MHz is the most common, we use only those data 
for this paper. Until 2009 January, the data were 
taken using the GBT spectrometer SPIGOT card 
with sampling time of 81.92 fis. SPIGOT provides 
1024 synthesized frequency channels across a 50 
MHz bandwidth. After 2009 January, the GBT 
spectrometer GUPPI was used. It covers 2048 
frequency channels with sampling time of 61.44 
/is, with a wider bandwidth of 200 MHz. All the 
data were dedispersed and folded using the pulsar 
analysis package SIGPROC, a ssuming a disper - 
sion measure of 48.914 cm ~ 3 pc (ILvne et al. 20041) . 



We used the ephemeris of ILvne et al.l (|2004l ) until 
2006, and sin c e then have used the ephemeris of 



Kra mer et al. I (|2006l) to form mean pulse profiles. 
A total of 52 data sets with 234 total hours of 
observations at 820 MHz have been analyzed. 

2.1. Pulse profile evolution of the two 
bright phases 

The pulsar is brightest in the orbital phase re- 
ion 185°-235°, and second brightest in 265° -305° 
Lvne et al.l liooi iBurgav et al.ll2005l ). Table 1 
shows the longitudes (as measured from the as- 



cending node) of each phase in which B is easily 
detected in early observations. In Fig. Q] and [21 
we present the integrated pulse profiles for BP1 
and BP2, respectively, on 12 days that pulsar B 
was observed at 820 MHz. The very first obser- 
vation, which was made in December 2003, shows 
a unimodal average pulse profile. This unimodal 
pulse profile gets broader and splits gradually into 
a two-component pulse profile over time. This 
evolution is common for both bright phases, but 
the specific features are somewhat different. 

The evolution from a single-peak profile to- 
wards a two-peak profile begins around 2006 May 
(MJD 53860). In BP1, the first peak is brighter 
than the second one. This two-peak profile ex- 
ists almost until the disappearance of the pulsar 
in this orbital phase region (see Fig. [T]). In BP2, 
the evolution has a similar pattern to BP1, but the 
remarkable difference is that the relative heights of 
the two peaks differ (see Fig. [2J . It initially looks 
similar to the profile of BP1, but on MJD 54055 
the relative amplitudes of the two peaks change 
and the second one becomes more prominent. This 
persists throughout the evolution and we believe 
that this feature is due to the impact of A's emis- 
sion on B, the effect of which changes over the 
period of the observation. Grey— scale plots of the 
first and the second bright phases of the two-peak 
profile are shown in Fig. [3] The first peak of each 
epoch has been normalized and aligned in both 
phases. It can be clearly seen that the two peaks 
are moving apart and the separation grows with 
time in both phases simultaneously. 

In order to get more information on the two 
peaks to describe the evolution, we fit two Gaus- 
sians separately for each peak of the two-peak 
pulse profile. From the Gaussian fits, we cal- 
culate the full-width-half-max (FWHM) of each 
peak of the two-peak pulse profiles in the two 
bright phases. They are shown in Figs. [4] and 
[SJ Linear least-squares fits qualitatively show that 
the FWHM of the first peak of BP1 decreases 
with time while that of the second peak increases 
slightly. In BP2, the FWHM of the first peak 
shows a slight decrease, but the decrease of sec- 
ond is more prominent. 

We present the evolution of the separation of 
the peaks in the two bright phases in Fig. |6lThe 
rates of separation are calculated from a lin- 
ear least-squares fit to be 2.6(1)° yr _1 and 
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Table 1: Starting and ending orbital longitudes of each phase in which B is easily detected. 





BP1 


BP2 


WP1 


WP2 


Start 


185° 


265° 


340° 


80° 


End 


235° 


305° 


30° 


130° 



Note. — The two bright phases are denoted as BP1 and BP2 and the two weak phases as WP1 and WP2. 



2.6(2)° yr" 1 for BP1 and BP 2, respectively. It 
appears that the profiles in both bright regions 
present the same rate of change in their compo- 
nent separation. 

2.2. Flux evolution of the two bright 
phases 

In both bright regions, the integrated pulse flux 
density has decreased gradually over time (see 
Fig. [T]and[2]). The pulsar was detected in both 
bright phases with the last significant detection in 
March 2008 (MJD 54552) at 820 MHz. We es- 
timated flux densities at 820 MHz using the ra- 
diometer equation. First we calculated the ra- 
diometer noise, using a system temperature of 
35 K (the system temperature is defined as the 
sum T sys = T roc + T spffl + T atm + T sky ), and con- 
sidered this as the flux density at the off-pulse 
region of the pulse profile. The flux density is 
then obtained by multiplying the pulse profile by 
the ratio of radiometer noise to the standard de- 
viation of the off-pulse phase and subtracting the 
mean off-pulse level. We have carried out this cal- 
culation for only the two bright phases, because 
the emission in the weak phases disappeared much 
earlier (discussed in §2.4 with more details). The 
calculated flux densities in both bright phases on 
MJD 52997 (0.95(4) and 0.65(4) mJy for BP1 and 
BP2 at 820 MHz, respectively) are con sistent with 
the va lue that has been calculated bv lLvne et al.l 
(|2004h (0-1.3(3) mJy at 1390 MHz). Fig.[7]shows 
the mean flux densities of different epochs which 
have been observed at 820 MHz. This confirms 
that the flux density significantly decreases over 
time and almost reaches zero around MJD 54852 
in both bright phases. The rate of decrease is cal- 
culated to be 0.177(6) and 0.089(7) mJy yr" 1 for 
BP1 and BP2, respectively. The flux densities of 
the last few epochs are only upper limits (denoted 



by arrows) and not included in the fits. Our tim- 
ing solution is not stable enough to provide a re- 
liable prediction of the expected phase of the pul- 
sar on these days, making it difficult to determine 
whether any apparent peaks are real. The peaks 
on MJDs 54856 and 54852 have the same pulse 
phase in both bright phases, suggesting they are 
real. However, given their low signal-to-noise, we 
describe them by upper limits. 

2.3. Analysis and comparison of light 
curves of the two bright phases 

The orbital-phase binned light curves of the two 
bright phases were obtained by integrating the flux 
in a window covering 5% of the spin period and 
centered on the pulse peak to reduce the effect of 
baseline noise. Then each light curve is smoothed 
by using a boxcar with a width of 30 pulses to 
reduce the significant pulse-to-pulse modulation. 
In Figs. [8] and [9l we present the light curves of the 
two bright phases on 12 different days that have 
been observed at frequency 820 MHz. Initially in 
both bright phases the emission was detected in a 
region covering ~ 30° of orbital phase, but this has 
shrunk over time, although this did not happen 
symmetrically in both phases. 

In order to analyze the evolution of the two 
bright phases with time, we calculated the start 
and end orbital longitudes, at the 10% of the 
maximum, of BP1 and BP2. These are shown 
in Fig. 1101 The linear least-squares fits show the 
start of BP1 moves to higher longitudes at a rate 
of 3.1(4)° yr _1 and the end of BP1 moves back- 
wards at a rate of -1.4(4)° yr _1 . In BP2, the both 
start and end orbital longitudes move to higher 
values at a rate of 4.1(2)° yr" 1 and 2.1(2)° yr" 1 . 
These rates show that BP1 and BP2 shrink at a 
rate of 4.5(6)° yr -1 and 2.1(3)° yr -1 , respectively. 
The above rates are updates for those presented by 
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Burgav et al.l ([20051 ) . These results confirmed that 
the disappearance did not happen symmetrically. 
However, the disappearance did occur simultane- 
ously in both bright phase regions. 

2.4. Analysis of the two weak phases 

Weak radio emission has been detected in two 
orbital phase ranges in addition to the two bright 
phases. One of the two weak emission regions, 
340° -30° (hereafter WP1) is much brighter than 
the other, 80°-130° (hereafter WP2), and is de- 
tectable until 2007 November (MJD 54429) (see 
Fig. EJ. WP2 vanished in 2006 February (MJD 
53783), more than a year before the disappearance 
of the WP1 (see Fig. [T2|) . In general, both weak 
phases show a unimodal pulse profile and no dra- 
matic pulse shape changes over time. Some data 
show two apparent peaks in their pulse profile, but 
the lower signal-to-noise makes it difficult to model 
these in the same way as for the bright phases (e.g. 
on MJD 53702 in WP1 and 53701 in WP2). If the 
evolution is solely due to relativistic spin preces- 
sion, we would expect the same profile evolution 
for all phases. However, the emission may not be 
strong enough in the weak phases to detect the 
second profile component. For instance, at MJD 
54113 the brightness of the second brightest com- 
ponent of BP1 is 30% of the brightest component 
and that of BP2 is 67%. In comparison, the root- 
mean-square noise in WP1 is almost 50% of the 
pulse maximum, making it difficult to detect the 
second peak of the profile. Therefore, within the 
uncertainties, the component separation evolution 
may be similar for bright and weak phase profiles. 

We calculated the mean flux densities of both 
weak phase regions (see Fig. [T3j) . In these phases, 
the mean flux densities are lower than that of 
bright phases. The decreasing rates of the mean 
flux densities are calculated to be 0.032(8) and 
0.02(2) mJy yr" 1 for WP1 and WP2, respectively. 
The larger uncertainty of the data points results 
in the poor quality of the fits. However, due to the 
weakness of the pulse emission in WP1 and WP2, 
we only use the two bright phase regions for our 
modeling. 

3. Determining the Geometry of the sys- 
tem 



Clifton fc Weisbergl (|2008h proposed a frame- 



work to model the 2-D pulse profile of a pro- 
cessing pulsar by using a simple circular hollow- 
cone beam. They applied this model for PSR 
B1913+16 to explain its long term pulse profile 
evolution. To represent the different intensity lev- 
els of the emission beam, they constructed a set of 
coaxial circular hollow cones with different angu- 
lar radii. Each circular cone represents a different 
intensity level. The outermost cone represents the 
lowest intensity level. The intensities gradually in- 
crease towards the maximum-intensity-level cone 
cone somewhere between the outermost cone and 
the center of the beam, and then decrease until 
reaching the center of the beam. With this special 
construction of the beam, it is possible to explain 
the hour-glass or oval shapes of the 2-D pulse pro- 
files. 

Moreover, this model constrains the geometri- 
cal parameters of the precessing pulsar. Mainly, 
the model predicts the longitudinal separation 
(i.e., pulse width) at a given epoch for a given 
intensity level for some set of geometrical param- 
eters. Then we fit this model-predicted width to 
the corresponding observed-pulse-profile width for 
a range of intensity levels to obtain the geometri- 
cal parameters. 

First we used the hollow-cone beam shape in 
order to determine the geometry of the B pulsar. 
Then we used a more complicated elliptical horse- 
shoe shaped beam, filled with different intensity 
levels as in the circular beam construction, to con- 
strain the geometry and found that this fit our 
observations better. 

In both cases, the required pulse-width data at 
different intensity levels of the pulse profiles have 
been calculated by fitting Gaussians to the profiles 
and then calculating the pulse widths at different 
intensity levels of the fitted Gaussian curves. The 
modeling has been done for the two bright phases 
separately to constrain two sets of geometrical pa- 
rameters. 

3.1. Modeling the geometry with a circu- 
lar hollow-cone beam 

In this section, we use a circular hollow-cone 
beam and the same set of equations as in CW08 
to determine the two angles a and 0, which are 
the magnetic inclination (i.e. angle between spin 
axis and magnetic dipole moment) and the colati- 
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tude of the spin axis (i.e. angle between spin axis 
and orbital angular momentum), respectively (see 
Fig. HH (a)). The model also includes another pa- 
rameter, To, the epoch for which the precession 
phase is equal to zero. The fit of the width at 
different intensity levels of our data to the model 
is shown in Fig. [T51 The fit has been done by 
minimizing the chi-square statistic over the en- 
tire parameter space of the two angles (a:[0,7r], 
(9: [0,tt]) and T : [1930,2001]. This range of dates 
encompasses one full precessional cycle for pul- 
sar B. The required errors of the data points for 
the chi-square statistic are obtained by the fitted 
Gaussians, which are the same for a given epoch 
at different intensity levels. 

The best-fit results for the two angles for BP1 
are a = 80.0°t{f Q ° and 9 = 20.0°t{f 5 o, while for 
BP2 a = 78.0°t^ 5 oo and 9 = 20.0°±i 2 °5 O . A sec- 
ondary peak in the PDFs caused large errors for 
the best fit parameters at their 68.5% confidence 
interval. The fits estimate T = 1959.4(1) yr for 
both phases. Thes e best fit values are very dif- 
ferent from those in lBreton et al.l (2008), who de- 
termined a = 70.9(4)° and 9 = 130.0(4)° (see Ta- 
ble 2). In addition to this, the simple circular 
beam shape predicts a disappearance in ~ 2018 
instead of 2008. Due to the poor fit and the in- 
ability to explain the disappearance of the radio 
emission in 2008, we experimented with other pos- 
sible beam shapes. 

3.2. Modeling the geometry with an ellip- 
tical horse-shoe shaped beam 

In this case, we used an elliptical hollow-cone 
filled beam with different intensity levels and then 
assumed that only a part of this beam is de- 
tectable, resulting in a horse-shoe shaped beam. 
To simplify the geometry, we assumed that the 
beam shape is fixed with respect to the neutron 
star and it rotates around the spin axis in a way 
such that the semi-major axis is always aligned in 
the radially outward direction from the spin axis 
(see Fig. [T6]) . In order to model this particular 
shape, we derived an equation for the profile width 



Wj(t) = arcsin[2 \J[n x {t)] 2 + [n y (t)] 2 smi]j(t)], 

(1) 

where Wj(t) is the width of the pulse profile at a 



given time, with the subscript j specifying the in- 
tensity level of the profile. The line-of-sight vector, 
h(t), can be expressed as in equation (3) of CW08 
and n x (t) and n y (t) in the above equation are the 
X and Y components of h(t). The third variable, 
r]j(t), is the angle which is subtended at the ori- 
gin of the coordinate system due to the encounter 
of the line-of-sight vector with the elliptical beam; 
namely, it is the angle between the line joining the 
origin with the point at which the line-of-sight vec- 
tor encounters the beam and the line joining the 
origin with the center of the elliptical beam (see 
Fig. [HI). Thus, amending the expression of CW08 
for an elliptical beam, rjj(t) can be expressed as 



r n 2 (t) 



. 2 

sin a 



[nit) 



sin a] 2 [1 - R 2 } 



2n(t) sin a 



(2) 

where Rj = bj/aj and a,j and bj are the semi- 
major and semi-minor axes of the elliptical beam 
for a given intensity level. They can take different 
values corresponding to the intensity levels of the 
beam in a way such that cij/bj, or the ellipticity 
of the beam, is a constant. 

The horse-shoe shaped beam is a subsection 
of the elliptical beam. We can construct this by 
restricting the detectable area of the beam to a 
section of the elliptical beam, with r% and r 2 the 
radial distances on the X-Y plane of the ellipti- 
cal beam within which emission is detectable (see 
Fig. HU (b)). Simply, we can detect the emis- 
sion beam of the pulsar only when sin a + n < 
\Z[n x (t)] 2 + [n y (t)] 2 < sin a + r 2 , where n x (t) and 
n v (t) are the X and Y components of n(t). 

The fit has been done for the two bright phases 
separately by searching over the entire parame- 
ter space for the two angles a and 9 and the 
ratio ajb using a like l ihood analysis (see, e.g., 
McLaughlin fc Corded (|2000h ). First the indi- 
vidual likelihood of each data point was calcu- 
lated by assuming a Gaussian distribution such 
that the likelihood for model 9 and measure- 
ment i is £i(8) = (27rcr, 2 ) _1/2 cxp(-(w modo i jl - 

Wmcasurocu) 2 / 2cr2 )> where Wmodcl.i and ^measured,* 

are the model and measured widths and <7i is the 
error on the measured width. We then multiplied 
all A (8) for all data points i (i.e. profiles from 
all epochs and at all intensity levels) to get the 
total likelihood £(0) for parameter combination 



b 2 



G 



0. Aot, the total likelihood over all parameter 
combinations, was calculated by summing all the 
likelihoods over all parameter combinations. Then 
the posterior probability distributions for param- 
eter values were calculated by dividing the likeli- 
hoods for a particular parameter value by the total 
likelihood £tot > assuming flat prior probability dis- 
tributions. The best fit 2-D pulse profile is shown 
in Fig. [17] and the posterior probabilities of the 
three parameters a, 9, and a/b for BP1 are shown 
in Fig. EES] 

The fitting constrained the two angles for BP1 
to be a = 65.0°ti$ and 9 = lte.5°±lfo, while 
for BP2 a = 73.0°t3;5l and 9 = 138.5°±*f a . We 
expect the same estimates for the two angles in 
both bright phases, indicating that our errors are 
underestimated or the model may not fit the data 
pe rfectly. These angles are similar to those derived 
bv lBreton et"aTI (|2008h (see Table 2). The best-fit 
estimate of the ratio of the beam semi-major and 
semi-minor axes (a/b) is 1.72^o'2o an( i ^-^-0.22 
for BP1 and BP2, respectively. Because of the 
varying influence of A, the magnetosphere could 
have different shapes at different orbital phases, 
and these angles could in fact differ. In order to 
cause the disappearance in 2008, the ratio r\jr% 
has to be 0.70 and 0.48 for BP1 and BP2, respec- 
tively. The parameter To is 1950.3(1) yr for both 
bright phases. Since the spin axis takes 75 yr to 
precess a full cycle around the orbital angular mo- 
mentum axis, the year 2008, within which disap- 
pearance occurred, is not a special year relative to 
this estimated To. 

4. Discussion 

Pulse profile changes have been detected in 
several othe r double neutron star systems: PSR 
B1913+16 (IWeisberg et all Il989i: IKramerl Il998l 
Weisberg fc Tavlorl2005l ). B1534+12 dArzoumanian 



1996; St airs et al.ll2004h J1141 -6545 (jHotan et al 



20051: iManchester et all l2010lh and J1906+0746 
( Lori mer et al J l2006h . In these pulsars, the long- 
term pulse profile evolution has been interpreted 
as the variation of the radio beam orientation 
with respec t to our l ine-of -sight due to geodetic 



precession. IKramerl (|1998l ) demonstrated a fully 



consistent model for PSR B1913+16 assuming a 
circular hollow-cone beam geometry and the pre- 
cession rate predicted by general relativity. More 



recently, CW08 revisited the problem with a dif- 
ferent geometrical framework using the circular 
beam approximation. 

We used the CW08 model in order to constrain 
the geometry of pulsar B. First we used the circu- 
lar hollow-cone beam shape they proposed in their 
paper and estimated the geometrical parameters a 
and 9 by fitting our data. The best fit values are 



different from the measurements of iBreton et al 



(|2008f ) and the fit is poor. It is likely that the dis- 
cr epancy betwe e n our best-fit geometry and that 
of IBreton et alJ ([2008), as well as the inability of 
the CW08 to fully explain our data, are due to a 
non-circular emission beam shape for pulsar B. 

In §3.2, we modified the geometrical frame- 
work of CW08 with an elliptical hollow-cone beam 
shape. This produced a much better fit than the 
circular beam. Furthermore, the two angles a and 
9 are consistent in the two bright phase s and with 



the predictions of IBreton et al.l (|2008f) (see Ta- 
ble 2). To obtain the disappearance in 2008, the 
detectable area of the beam has to be restricted to 
a section of the elliptical beam, producing a horse- 
shoe shaped beam. Given this model, we calculate 
the ratio of the semi-major and semi-minor axes 
(a/b) to be 1.72±°;ig and l.lbt° Q f 2 for BP1 and 
BP2, respectively. After transforming this ratio 
to the assumed unit sphere by dividing by cos(a), 
we estimate the ellipticity of the radio beam of B 
in the two bright phases separately: 0.9(1) and 
0.9(2) at BP1 and BP2, respectively. The ratio 
r\/ri in the two bright phases reveals that only 
29(6)% of the beam area was detectable in BP1 
and only 21(5)% in BP2. The different estimates 
for the effective area of the beam in the two bright 
phases show that the emission beam may have dif- 
ferent orientations corresponding to the different 
orbital phases. The characteristic pulse shapes are 
also different at different orbital phases. This is a 
direct consequence of the influence of pulsar A's 
wind on the em ission beam of pulsar B (see, e.g., 
Lvutikovl (|2004j )V 



The special choice of the beam construction 
with different intensity levels explains the single- 
peaked to double-peaked pulse profile evolution. 
On the first day, MJD 52997, the observed bright 
single-peaked pulse profile can be explained as 
our line-of-sight just grazing the outer edge of 
the maximum intensity level of the beam. When 
the line-of-sight gradually moves inward with time 
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Table 2: Geometrical parameters of pulsar B deriv ed from our circular beam model (CBM) and elliptical 
beam model (EBM) and the eclipse model fitting by Breton et al. ( 2008f) . 









O 




en 


CBM 


:BP1 


80.0(- 


-15.0, +1.5) 


20.0(- 


16.5, +1.5) 


CMB 


:BP2 


78.0(- 


-11.0, +1.5) 


20.0(- 


12. 5, +1.0) 


EBM 


:BP1 


65.0(- 


-1.0, +2.0) 


134.5( 


-6.0,+1.5) 


EBM 


:BP2 


73.0(- 


-3.5, +1.4) 


138. 5( 


-7.5,+4.0) 


Breton ct al. (2008) 


70.9(- 


-0.4, +0.4) 


130.0( 


-0.4,+0.4) 



(i.e., towards the spin axis; see Fig rj+T b)). the 
two-peaked pulse profile is formed since the max- 
imum intensity level of the beam crosses our line- 
of-sight twice. However, the observed decreasing 
flux density over time cannot be explained by ei- 
ther model, since the intensity is constant along 
the cross-section of a hollow cone of the conical 
beam, in both the circular or elliptical cases. We 
believe the observed decrease in flux density must 
be due to a gradient in brightness across the beam 
or due to the changing influence of A on B dur- 
ing the span of these data. The distance between 
the two pulsars changes over time relative to our 
line-of-sight due to the relativistic advance of pe- 
riastron of 17° yr" 1 of the system. This results in 
a different amount of impact on B's emission by A 
with respect to Earth and leads to a change in the 
shape of the horse-shoe beam slightly, but this is 
likely a small effect, and not easy to model. Thus 
we ignored this effect. 

The data show that the pulse profile had two 
peaks when the disappearance occurred in 2008 
and clearly did not become a single peak before 
the disappearance. The horse-shoe shaped beam 
model can explain this phenomenon qualitatively, 
as the line-of-sight moved out of the elongated side 
of the beam (i.e. the two opened arms of the horse- 
shoe beam, see Fig. [14] (b)). This results in the 
disappearance of both peaks of the pulse profiles 
in both bright phases simultaneously. The reap- 
pearance of the emission of pulsar B is predicted to 
occur in ~ 2035, according to the horse-shoe beam 
model, assuming emission from the same part of 
the beam. This is because our line-of-sight en- 
counters the emission beam twice during the time 
that the spin axis rotates around the orbital an- 
gular momentum axis once (~71 yr). This can 
be understood by visualizing the geometry of the 



emission beam with two angles, a ~ 65° and 9 « 
130°. Thus the spin axis of the pulsar B is below 
the inclination plane and the magnetic moment 
axis crosses the inclination plane twice during ev- 
ery 2.8 sec period. Since the inclination plane of 
the system is almost edge-on with respect to the 
Earth, we are able to detect the emission beam 
twice during the 71 yr period. Geodetic preces- 
sion along with the unique geometry of the sys- 
tem is responsible for this phenomenon. Note that 
this prediction for the reappearance time is only 
valid assuming a single horse-shoe emission com- 
ponent. If the beam is symmetric, with two horse- 
shoe shaped emission components separated by a 
gap, then we would expect a reappearance around 
2014. Because we have only sampled a limited 
part of the beam, it is impossible to know what 
the true intensity distribution of B's beam is at 
this point. 

Despite the unusual modulation of B's emis- 
sion by A, the 3-5% duty cycle of B's pulse and 
the pulse profile shapes throughout the evolution 
are similar to those of pulsars with similar peri- 
ods, with the dramatic geodetic precession of B 
offering us the chance to sample a much larger 
portion of the beam than typically possible. It 
is difficult to determine what is responsible for 
B's flux density variation. It may be that most 
pulsars have instrinisic intensity variations across 
their beams. Or, the decrease in the intensity of 
B's emission with time may be solely due to B's 
magnetosphere shape changing with the varying 
influence of A. Theoretical modeling of the influ- 
ence of A's wind with time is necessary to address 
this. Our horse-shoe shape d beam f its th e stan- 
dard conal emission model iRankinl (|1983I ). with 
the observed single to double peak evolution due 
to sampling different lines-of-sight across the con- 



S 



ical beam, but this model does not explain the 
decreasing flux density and eventual disappear- 
ance of the pulsar. Pulse profile shapes have 
also be en explained through the 'patchy' beam 
model ( Lvne fe; Manchester! 19881) in which inten- 
sities across the beam vary widely and often only 
a portion of the cone is visible. Observations 
of relativistic binary pu lsar system J 1141— 6545 



( Manchester et al. 2010l ) show that the emission 



across the beam is asymmetric with respect to the 
magnetic axis, consistent with a patchy beam. It 
is more difficult to interpret the evolution of B's 
profile with the patchy beam model, however, un- 
less the patches are large enough to explain both 
components disappearing at the same time. Con- 
tinued monitoring and sensitive searches for emis- 
sion from B will provide invaluable information on 
the true beam shape of this pulsar. 

Lvutikovl (|2005h proposed a model to explain 



the change in the orbital longitudes of the two 
bright phases. It relies on the wind of pulsar A 
affecting the direction of emission of pulsar B and 
making it miss our line of sight for most of the 
orbit. In this model, the two bright phase re- 
gions evolve as a decreasing or increasing orbital 
phase gap between the two regions according to 
the direction o f the spin precession (see Fig. 6 of 
Lvutikovl 120051 ). However, the light curves of the 
two phases show that the real situation is some- 
what different (see Fig. fTOjl . Motion of the ending 
orbital phase of BP1 to lower phases and start- 
ing orbital phase of BP2 to higher phases result in 
an increase in the orbital phase gap between the 
bright regions instead of the motion of both the 
starting and the ending phases that is predicted 
by the model. We are working on a more quanti- 
tative model in order to explain the evolution and 
the disappearance of B's emission. 

We thank Joel Weisberg, Duncan Lorimer, 
and Maxim Lyutikov for very useful discussions. 
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Fig. 1. — The pulse profiles of BP1 on 12 different days including the very first observation (MJD 52997). 
All data have been observed at a frequency of 820 MHz. Each profile covers 20 min of orbital longitude (from 
185°— 235°) and there are 256 bins across the entire profile. Since predictions of absolute pulse phase are not 
available for these observations, the highest profile peak has been aligned to phase 0.5 at each epoch. The 
horizontal solid and dotted lines show the baseline, or off-pulse mean, of the profile and the corresponding 
standard deviation of the off-peak region, respectively. The effective time resolution of the profiles is 0.01 
sec. 
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Fig. 2. — As described in Fig. [Q but for BP2. 





Fig. 3. — Grey— scale plots of BP1 (left) and BP2 (right). There are 256 bins across each pulse profile. The 
first peak of each pulse profile has been normalized to unity and aligned to visualize how the second peak 
evolves. The discontinuities of the plots are artifacts due to boundary level changes. 
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Fig. 4.— Full width at half max (FWHM) of the first peak of BP1 (left) and the second peak of BP1 
(right). The linear least-squares fits show a decreasing FWHM rate of 0.59(8)° yr^ 1 for the first peak and 
increasing FWHM rate of 0.4(3)° yr" 1 for the second peak. 
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Fig. 5.— FWHM of the first peak of BP2 (left) and that of second peak of BP2 (right). The linear least- 
squares fits show a decreasing FWHM rate of 0.8(2)° yr _1 and 1.7(2)° yr _1 for the first and second peak, 
respectively. Note that due to the small number of bright pulses included in these composite profiles, the 
profiles show a high amount of variability that does not follow any trend. 




Fig. 6. — The varying separation of the two peaks of pulse profiles of BP1 (left) and BP2 (right). The linear 
least-squares fits show a separation rate of 2.6(1)° yr -1 and 2.6(2)° yr _1 for BP1 and BP2, respectively. 
Note that in BP1, the two-peak pulse profile lasts longer than in BP2. 
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Fig. 7. — Mean flux density of the radio emission in BP1 (left) and in BP2 (right). The data marked by 
the arrows are the upper limits and ignored in our linear least-squares fit. The decreasing rate of the mean 
flux density is calculated to be 0.177(6) and 0.089(7) mJy yr" 1 for BP1 and BP2, respectively. 
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Fig. 8. — The light curves of BP1 on 12 different days at an observing frequency of 820 MHz. Each data set 
has been cleaned to get rid of radio frequency interference. All light curves have been smoothed by using a 
boxcar averaging technique with a width of 30 pulses. The horizontal solid and dotted lines show the baseline 
of the plot and the corresponding standard deviation of the off-peak region (140° — 180° and 240°— 260°), 
respectively. The last data set shows just noise. 



14 




260 270 280 290 300 310 260 270 280 290 300 310 260 270 280 290 300 310 260 270 280 290 300 310 
Orbital phase Orbital phase Orbital phase Orbital phase 



Fig. 9.— As in Fig. HI but for BP2. 




Fig. 10. — Evolution of the start and end orbital longitudes of BP1 and BP2, measured at the 10% of the 
maximum. 
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Fig. 11. — As in Fig. [TJ but for WP1. Each profile covers 20 min of orbital longitude (from 340°-30°) and 
there are 256 bins across the entire profile. 
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Fig. 13. — Mean flux density of the radio emission in WP1 (left) and in WP2 (right). The data marked by 
the arrows are the upper limits and ignored in our linear least-squares fits. The decreasing rate of the mean 
flux density is calculated to be 0.032(8) and 0.02(2) mJy yr" 1 for WP1 and WP2, respectively. 
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Fig. 14. — Geometry of pulsar B with its spin and the precession in a Cartesian coordinate system with a 
unit sphere at the center; all of the trajectories move on this sphere, (a) The coordinate system is chosen 
with the Z-axis parallel to the spin-axis of the pulsar (£l sp in) with the X-axis in the plane of the spin-axis and 
the orbital angular momentum axis (fl or bit)', i-e., the & S pin and the £l rbit are always in the X-Z plane. The 
pulsar-earth line-of-sight is h(t) and dashes represent the conical trajectory of it on the unit sphere about 
&orbit due to geodetic precession. Dots represent the conical trajectory of the co-rotating magnetic field (ft) 
on the unit sphere about il sp in due to the pulsar spin. The inclination of the orbit is i. The colatitude 
of spin axis and the misalignment of the magnetic field are 9 and a, respectively, (b) The projections of 
trajectories of ft and h(t) on the X-Y plane. Dash-dot lines represent the different intensity levels of the 
elliptical beam. The intensity increases gradually from the outer most edge of the beam to inwards until the 
maximum intensity level (the solid line of the beam), and then decreases towards the center of the beam. 
The lightly shaded region represents the detectable beam area when the beam sweeps about the spin axis. 
This area is determined from the radial lengths r\ and r2, and they are measured from the center of the beam 
on the X-Y plane. With this restriction, the effective beam has a horse-shoe shape and it is represented as 
the dark shaded region. Dashes show the trajectory of the h(t) on the X-Y plane. When h(t) crosses the 
beam area (lightly shaded area), we are able to detect the pulsar. 
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Fig. 15. — The hour-glass shaped two-dimensional pulse profile of PSR J0737— 3039B, assuming a circular 
hollow-cone beam. The best fit parameters for BP1 are a = 80. ^^ and 6 = 20.0 o t\f 5 o (left), and for 
BP2 arc a = 78.0°i 11 5 o and 8 = 20.0 o ti 2 °5 o (right) . Fitting has been done by searching the entire range 
of the two angles. The dots are the widths at equal intensity of the pulse profile. Each horizontal row of 
dots represents an observation at a given epoch. Equal-intensity contours are produced from the emission of 
the circular symmetric conical beam. The solid lines are the equal intensity level of the peaks of two-peak 
pulse profile. The intensity increases from the inner dashed line outwards until the first solid line, which is 
the intensity of the peak, and then decreases outwards again. The intensity levels are, from inner dashed 
line, 80%, 90%, 100%, 90%, 80%, 70%, 60% and 50%. The vertical axis is calibrated in years, and can also 
be considered the spin precession phase, which increases linearly with time. The time of precessional phase 
zero in the fits is T = 1959.41 yr for both phases. We have omitted the data between MJDs 53175 and 
53939. During this time, the profile was bi-modal but the separation between the peaks was small, making 
it difficult to measure widths at different intensity levels. 
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Fig. 16. — A zoomed-in view of the elliptical beam on the X-Y plane with a constant intensity level. The 
given constant intensity level crosses our line-of-sight when the line-of-sight is within the effective beam area 
(lightly shaded area in Fig I14f b)). The semi- major and semi- minor axes, Uj and bj, can take any value 
corresponding to the intensity level in a way such that the ratio &j/hj is a constant. After projecting to the 
X-Y plane, the semi-major axis is aligned in the radially outward direction from the spin axis, sin(a) is the 
distance to the center of the beam from the spin axis on the X-Y plane, rjj(t) is defined also on the X-Y 
plane, and it is the angle between the line joining the origin with the point at which the line-of-sight (h(t)) 
encounters the constant intensity level of the beam and the line joining the spin axis with the center of the 
elliptical beam. 
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Fig. 17. — The two-dimensional pulse profile of PSR J0737— 3039B, assuming a horse-shoe shaped beam. 
The best fit parameters for BP1 are a = 65.0 o j^;g°, 6 = 134.5°t^°, and a/b = 1.72±g;^g (left). For BP2, 
the best fit parameters are a = 73. !^;^, 9 = 138. 5°^f®o, and a/b = 1.15^0.22 (right) . We have omitted 
the data between MJDs 53175 and 53939. During this time, the profile was bi-modal but the separation 
between the peaks was small, making it difficult to measure widths at different intensity levels. 




Fig. 18. — The calculated posterior probabilities of the parameter values given by likelihood analysis for 
BP1. The best fit parameters and their 68% confidence intervals are estimated to be a/b = 1.72+Q20i a = 
65.0°l?;go and = 134.5°±J;&. 
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